Inelastic Scattering from Core-electrons: a Multiple Scattering Approach 
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The real-space multiple-scattering (RSMS) approach is applied to model non-resonant inelastic 
scattering from deep core electron levels over a broad energy spectrum. This approach is applicable 
to aperiodic or periodic systems alike and incorporates ab initio, self-consistent electronic structure 
and final state effects. The approach generalizes to finite momentum transfer a method used exten- 
sively to model x-ray absorption spectra (XAS), and includes both near edge spectra and extended 
fine structure. The calculations can be used to analyze experimental results of inelastic scattering 
from core-electrons using either x-ray photons (NRIXS) or electrons (EELS). In the low momentum 
transfer region (the dipole limit), these inelastic loss spectra are proportional to those from XAS. 
Thus their analysis can provide similar information about the electronic and structural properties of 
a system. Results for finite momentum transfer yield additional information concerning monopole, 
quadrupole, and higher couplings. Our results are compared both with experiment and with other 
theoretical calculations. 

PACS numbers: 79.20.Uv, 78.70.Ck, 71.15.Qc 
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I. INTRODUCTION 



The primary aim of this work is to model inelastic scat- 
tering quantitatively using a generalization of the real 
space multiple scattering (RSMS) approach. This ab ini- 
tio Green's function method has been extensively devel- 
oped for calculations of x-ray absorption spectra (XAS) 
and related spectroscopies Q], Q- Thus the approach 
permits calculations of inelastic scattering in arbitrary 
aperiodic materials over a broad energy range, including 
the near edge and extended fine structure. In inelas- 
tic scattering experiments, the scattering cross section is 
measured as a function of the energy to and momentum 
q transfered from the probe to the system (throughout 
this paper we use Hartrec atomic units, e = K = m = 1). 
When the energy transfer is close to the binding energy 
of a core state, the inelastic scattering cross-section ex- 
hibits pronounced jumps, from which one can study the 
core-excited states of the system. We focus in this paper 
on the excitations from such core levels. 

In the low momentum transfer regime the scattering 
cross section is dominated by the dipole allowed exci- 
tations. Qualitatively the dipole approximation is valid 
when qa <C 1 , where a ~ l/£ is the mean radius of the 
core state, ( being the effective core charge. In the dipole 
limit the inelastic scattering cross section is proportional 
to the x-ray absorption coefficient 0, 0, , and the di- 
rection q of the momentum transfer vector plays the role 
of the XAS polarization vector e. Thus the anisotropy of 
the x-ray edge can be studied by varying q in much the 
same way it is studied in XAS by varying the direction 
of e. However, as the magnitude of q increases, contribu- 
tions from other (dipole forbidden) excitation channels 
become important. The symmetries of a system gener- 
ally restrict the available excitation channels at a given 
excitation energy. The relative weights of the various ex- 
citation channels, as a function of momentum transfer, 



also depend on the spatial extent of the excited state. 
In general the classification of the spectral features into 
given spatial symmetries of the final state electron (i.e. 
s-type, p-type etc.) is not as straightforward as in XAS. 
This is because contributions to the spectra from different 
excitations often overlap at finite momentum transfer. 
Due to this overlap, a detailed analysis of experimental 
results can be difficult in the absence of quantitative the- 
oretical calculations. There are some exceptions though; 
core-excitons can have well defined spatial symmetry and 
energies well separated from other excitations. 

Non-resonant inelastic scattering spectra can be mea- 
sured using either electron energy-loss spectroscopy 
(EELS) jol or non-resonant inelastic x-ray scattering 
(NRIXS) 0- The cross sections in both cases are re- 
lated to the dynamic structure factor 5(q, uj) of a system. 
For the case of x-rays, the non-relativistic Born approx- 
imation for the NRIXS double differential cross-section 
is HQ 



(1) 



Here (da / d£Y)Th is the Thomson scattering cross-section. 
The Thomson cross section can be expressed in terms of 
the incoming (scattered) photon polarization vector and 
energy {ei,u>i} ({£ 2 ,^ 2 }) as 



d<r\ 2 „ 2 u> 2 

-jo =rS(ei ■e 2 y — , 
d\l J Th oji 



where r$ is the classical electron radius. In Eq. Q wc 
have implicitly assumed that the energy of the incident 
x-ray is not close to any core binding energy JSj (typically 
JSj <CUi). This is in contrast with resonant inelastic x- 
ray scattering experiments where the incident x-ray en- 
ergy is tuned close to a core binding energy and resonant 
processes dominate the double differential cross section. 
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The different forms of inelastic x-ray scattering and their 
relationships are discussed in more detail in Refs. H, @- 
In EELS the corresponding scattering cross section in the 
Born approximation is [lOj 

d^ = {dn) e J M > ^ 

where {da / <Kl) e - e = q~ i (p' jp) is the electron-electron 
scattering cross-section with p (j>') being the momentum 
of the incoming (scattered) electron. In other words, with 
these approximations the cross-sections differ only by a 
coupling constant which is related to the probe (i.e., a 
photon or electron) that interacts with the system. The 
system dependent part, namely the dynamic structure 
factor S'(q, to), is the same for both NRIXS and EELS, 
and depends only on the structure and excited state prop- 
erties of the system under investigation. 

Non-resonant inelastic x-ray scattering from core- 
excited states is often called x-ray Raman scattering 
(XRS). Traditionally XRS experiments have been limited 
to if -edges of low- .2? materials. However, for the if-edges, 
experimentally accessible systems include such important 
cases as those with carbon and oxygen, as has already 
been demonstrated (see Ref. ^2 an d references therein). 
Recently XRS experiments were done for the L-edges and 
./V-edges of materials with higher Z including N-edge of 
Ba in BasSi46 [TJ. The smallncss of the x-ray scattering 
cross-section implies that bulk information is generally 
obtained in these probes, i.e., that multiple scattering of 
the probe particles and surface effects are usually not se- 
rious problems. XRS has been used to study both the 
momentum transfer magnitude fLU Ibil ITU \wL Il7j and 
direction dependence at if -edges. Recently the K- 
edge of Be metal was studied using both degrees of free- 
dom (magnitude and direction) 18]. The typical energy 
resolution of current XRS experiments is around 1.0 eV 
FWHM, although some experiments with a energy res- 
olution of 0.4 cV have been carried out The finite 
core-hole lifetime broadening also limits the useful energy 
resolution to a few tenths of an eV. 

As noted above, one can also measure the fine struc- 
ture of the spectra above deep core edges using electrons 
as a probe. Unlike non-resonant x-ray scattering, EELS 
is not limited to \ow-Z materials, and has a typical en- 
ergy resolution on the order of 0.1 eV. Electrons interact 
strongly with solids, so that one is often able to obtain 
results with good statistical accuracy. However, their 
strong interaction also implies that multiple scattering 
of the probe electrons can be a problem which has to 
be removed or otherwise accounted for in the analysis 
of experimental results. This is especially true of ex- 
periments with large scattering angles or low incident 
electron energies, in which case one should explicitly in- 
clude m ultip le scattering in the analysis, as explained in 
Ref. [HHf. The dependence of x-ray edges on the mag- 
nitude of the momentum transfer has been investigated 
in Ref. HElil, and more recently in Ref. ■ The de- 
pendence on the momentum transfer direction at x-ray 



edges has also been used in the analysis of anisotropy 
of x-ray edges (see Ref. and references therein). 

Additionally the very high spatial resolution available in 
EELS makes it possible to study the anisotropy of x-ray 
edges in fine-grained solids |25j . 

Theoretically the problem of calculating inelastic scat- 
tering spectra from core-excited states is analogous to 
the calculation of core level x-ray absorption, the main 
difference being in the nature of the transition matrix 
elements. Thus there are number of effects that must 
be taken into account for accurate computational re- 
sults. Within the independent electron approximation, 
the main problem is the behavior of the core-excited fi- 
nal states. This can be approached cither by calculating 
the core-excited states one by one (e.g., in band-structure 
or molecular-orbital based approaches) or equivalcntly, 
by calculating the one-electron Green's function for the 
excited states, as is done in this work. For accurate re- 
sults one needs a highly quantitative method for calcu- 
lating single-particle properties such as electron densities 
and Coulomb potentials. Also a number of many-body 
properties of the excited states of the system have to 
be considered. These include the interaction of the final 
state electron with the core-hole created in the excita- 
tion process, as well as an approximation for the energy- 
dependent self-energies H(E) of both the electron and 
the core-hole, e.g. to account for inelastic losses. A 
good approximation in many cases is to calculate the 
final state in the presence of the core hole and the initial 
state with ground state potentials This approxima- 
tion is referred to as the final state rule. In principle 
one should include the whole energy dependence of the 
self-energy, as discussed for example in Ref. [2(|. How- 
ever, often a quasiparticle approximation is used both 
for the electron and the core-hole. In the quasiparti- 
cle approximation, the imaginary part of the self-energy 
causes the electron (and the hole) in the excited state 
to have a finite lifetime. Additionally the real-part of 
the self-energy shifts the spectral features from the posi- 
tions predicted by ground-state (or mean-field) approxi- 
mations, whereas traditional ground state density func- 
tional treatments can lead to significant errors in peak 
positions and intensities. For a full treatment, especially 
at high energies, one also needs to consider coupling to 
vibrational degrees of freedom, which give the spectra a 
temperature dependent damping comparable to that of x- 
ray Dcbye- Waller factors. A more detailed review of the 
problem of computing core-excited states in XAS can be 
found in Ref. [l| . As noted above the main difference be- 
tween inelastic scattering and x-ray absorption lies in the 
transition matrix elements. This is due to the fact that 
the operator mediating the transitions in inelastic scat- 
tering is cxp(iq ■ r), as compared to the dipole operator 
e-r (and, in some cases, also the quadrupolc operator) in 
x-ray absorption. We will discuss this difference in more 
detail below. 

In this article we have applied the real-space multiple- 
scattering (RSMS) approach to calculate the core-excited 
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states in EELS and NRIXS over a wide energy range. 
The RSMS approach adopted here, is essentially an ef- 
fective, independent particle approach for excited states 
which takes into account the most important final state 
effects. The approach is an extension of the RSMS ap- 
proach previously used for x-ray absorption spectra 
The momentum transfer dependence of both the near 
edge spectra and the fine structure in inelastic loss spec- 
tra are analyzed in detail and compared with experimen- 
tal results. 

Recently a two-particle approach to core-excited states 
has been developed [28ll29l | based on the Bethe-Salpeter 
equation (BSE). The BSE goes beyond the independent 
particle approximation with an explicit treatment of the 
screened electron core-hole interaction. However, when 
the core-hole interaction is stron g, th e results from the fi- 
nal state rule can be comparable [27]]. The BSE approach 
has been applied to different situations using the momen- 
tum transfer dependent matrix elements in the analysis 
of XRS HI III 13 For completeness our calcula- 
tions are also compared with calculations based on the 
BSE. However, the BSE approach is currently applicable 
only relatively close to a given edge, typically between 
a few tens and 100 cV depending on the material, and 
becomes computationally impractical for treating the ex- 
tended fine structure in inelastic scattering spectra. 

II. REAL SPACE MULTIPLE-SCATTERING 
FORMALISM 

A. Green's function formulation 

We now briefly describe the extension of the real- 
space multiple-scattering approach for calculations of 
non-resonant inelastic scattering. The RSMS approach 
can be regarded as the real space analog of the KKR 
band structure method [3(| , but unlike KKR, the method 
makes no assumption of symmetry and is applicable to 
periodic and aperiodic systems alike. Moreover, the im- 
plementation used here is an all-electron approach which 
can be applied to arbitrary systems throughout the peri- 
odic table. This method has been extensively developed 
for ab initio calculations of x-ray spectra including both 
EXAFS (extended x-ray absorption fine structure) and 
XANES (x-ray absorption near edge structure) 0, ■ Its 
generalization here to finite momentum transfer q is rel- 
atively straightforward, though non-trivial in several re- 
spects, as discussed below. The contributions to 5(q, uj) 
from excitations from a tightly bound core state \i) can 
be approximated using Fermi's Golden rule, i.e., 

Sfaw) = Yl l(/|e 4q - r N)| 2 <5(^ + Ei- E f ), (3) 
/ 

where Ei (Ef) is the initial (final) state quasiparticle en- 
ergy of the electron. For the remainder of the paper we 
will drop the index i and use simply 5(q, uj) when refer- 
ring to the contribution from a core state \i). Within the 



one-electron approximation and the final state rule, the 
final (photo-electron) states |/) are quasi-particle states 
which are eigenfunctions of the final state Hamiltonian in 
the presence of an adiabatically screened core hole (de- 
noted with a prime), i.e., H' = p 2 /2m + F c ' oul + T,(E). 
Here £(-E) is the self-energy (or dynamically screened 
exchange-correlation potential) which accounts for in- 
clastic losses, which are essential for a proper treatment 
of inelastic electron scattering. In our calculations we 
use the energy dependent local density approximation for 
£(£■) of Hedin and Lundqvist [3l| , based on the plasmon- 
pole approximation for the dielectric function. This ap- 
proximation is usually adequate at moderate to high pho- 
toelectron energies [lj. The states \i) are deep-core lev- 
els of the initial state Hamiltonian without a core-hole. 
For small momentum transfers in the case of core-excited 
states, the excitation operator exp(iq-r) can be expanded 
as 

exp(iq • r) w 1 + iq • r + 0(q 2 ). (4) 

Thus the dipole selection rule is approximately valid at 
small q, since the first term should not contribute to tran- 
sitions due to orthogonality of the initial and final states, 
and the next term (cx q 2 ) is also small. In fact since ui 
and q can, in principle, be chosen separately, the dipole 
approximation can sometimes be satisfied better in in- 
clastic scattering than in absorption, where the photon 
momentum and energy are always linked. At higher q, 
dipole forbidden transitions become important, and it is 
more convenient to expand the exponential in terms of 
spherical harmonics 

cxp(iq-r) = 4nJ2i l j l (qr)Y l * m (q)Y lm (r). (5) 

In) 

Thus for any finite momentum transfer all excitation 
channels are present . 

To avoid explicit calculations of final states it is ad- 
vantageous to re-express the Golden rule in terms of 
the one-particle Green's function or propagator G = 
[E — H' + iT}^ 1 in real space, where T is the core- 
hole lifetime. Thus using the spectral representation, 
-(l/7r)ImG(£0 = X f \f)(f\5(E - E f ) = p(E) the dy- 
namic structure factor becomes 

S(q, u) = (i\e-^ r 'Pp(r>, r, E)Pe^ r \i), (6) 

where the photoelectron energy is E = uj + Ei. The 
operator P projects the Green's function on to the unoc- 
cupied states of the initial state (without the core hole) 
Hamiltonian j2||. The operator P is needed here since, 
in general, the eigenstates of the final state Hamiltonian 
are not strictly orthogonal to the initial state of the scat- 
tering. In practice this non-orthogonality implies that at 
low momentum transfer 

(i\e- i ^\f)^(i\f)-i(i\c l -v\f) + ---. (7) 

Thus the dipole forbidden excitations, which incorrectly 
vary as oc (i\f), can begin to dominate over the dipole 
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allowed transitions oc q as q — > 0. The application of 
the full projection operator P would require taking an 
inner-product and appropriate matrix elements of all the 
occupied states \j) of the initial state Hamiltonian, i.e., 

occ 

Pe i **\i)=e i **\i)-^2\j)(j\e i * r \i). (8) 

3 

Instead we simply approximate the effect of P by modi- 
fying the excitation operator as 

Pe iq ' r |£) « [e iq ' r - {i\e l * r \i)]\i). (9) 

Clearly at low momentum transfers this gives matrix el- 
ements 

(i|e-^ r P|/)«-i(i|q-r|/)+e(g 2 ), (10) 

while at high momentum transfers the correction 
(i\e zcl ' r \i)(f\i) approaches zero. Although in what follows 
we do not explicitly include the operator P, it is included 
in the actual calculations. The terms coming from the 
difference of P and P are neglected. Although this is a 
somewhat uncontrolled approximation, comparison with 
experiments and previous XRS calculations at different 
momentum transfers show that these terms should be 
small, as illustrated in Fig. [21 below. Also work in XAS 
suggests that the difference not significant compared with 
other theoretical uncertainties [26j . 

In multiple scattering (MS) theory, the scattering per- 
turbation is the total electron potential, which is sepa- 
rated into contributions from "scattering potentials" vr 
which are localized on each atomic site R, i.e., 

K' oul + £(£) = ]>>A(r-R). (11) 

R 

In the RSMS method vr(t) is usually taken to be spher- 
ically symmetric. This is a good approximation for elec- 
tron scattering calculations at moderate electron kinetic 
energies, i.e., a few eV above the threshold. These po- 
tentials are calculated self-consistently by iterating the 
total electron density, potential and Fermi-energy, typi- 
cally requiring about 10-20 iterations. Once the poten- 
tials are known, scattering phase shifts 5i are calculated 
and dimensionlcss i-matrices evaluated using the relation 
ti = cxp(i<5;) sin((5/). With spherically symmetric poten- 
tials, the propagator G(E) (and hence the density ma- 
trix p(E)) can be represented in an angular-momentum 
L = (i, m) and site R basis |L,R). Thus at the site of 
the core-excited atom (R = 0) 

p(r', r 1 E) = Y, Rl(?, E)p LtL ,{E)R L ,{v' , E), (12) 

L,L' 

where Rl{y,E) are scattering states at energy E, which 
are regular at the origin. The expansion for r and r' 
about different sites R and R' is similar, with Gl,l'(E) 
replaced by Glr,l'R'(E). Consequently the calculation 



of the dynamic structure factor is reduced to a calcu- 
lation of (embedded) atomic transition matrix elements 
M L (q,E) = {R L (E)\e ici - r \i) and a multiple-scattering 
matrix p L ,L>(E) = Plo,l'o(E), i.e., H2 

S(q,w) =J2 M L(-<l,E)p L , L ,(E)M L ,(q,E). (13) 

LL' 

Here pl,u{E) = (—1/tt)IuiGl,l'(E) denote matrix el- 
ements of the final state density matrix, including the 
effect of the core-hole potential. It can be shown that 
5(q, ui) also satisfies a generalized oscillator strength 
sum-rule. The representation in Eq. (|13fl shows that the 
essential physics of the problem separates into two parts: 
i) a q dependent transition matrix which governs the pro- 
duction of photoelectrons into various final states, and ii) 
a propagator matrix Gl,l' which describes the scattering 
of the photoelectron within the system at a given exci- 
tation energy. The transition matrix elements Ml (q, E) 
are calculated using the expansion in Eq. (5) of exp(iq-r) 
in terms of spherical harmonics. Depending on the ex- 
citation energy and momentum transfer, different terms 
L in this expansion are important. For low momentum 
transfer and excitation energy the small Al transitions 
are most important, starting from the dipole Al = ±1, 
monopolc Al = 0, and quadrupole Al = {±2, 0} tran- 
sitions. Thus provided one can neglect the coupling to 
phonons (which is a good approximation for core exci- 
tations discussed in this work), a typical spectrum at 
low q can be analyzed using only these three excitation 
channels. Monopole transitions, which are present in 
S(q,u) through the term cx 4irj (qr)Yo Q (q)Yoo(r), have 
no counter-part in absorption (s-to-s transitions being 
forbidden within the dipole approximation). When the 
momentum transfer is increased in inelastic scattering, 
other excitation channels become more important, espe- 
cially at high energy values. 

The above RSMS formulation can be advantageous 
even for crystals, since periodicity is broken by the core- 
hole interaction and spectral broadening from the core- 
hole lifetime and the self energy. Typically Y c h + |Im£| 
is several eV at photoelectron energies above about 30 
cV of an edge. This broadening limits the range probed 
by the photoelectron to clusters of order a few hundred 
atoms, and gives rise to a short range order theory for 
energies above about 30 eV. Thus long range effects such 
as the sharp van Hove singularities of band structure cal- 
culations are naturally smeared out. Moreover, in this 
extended energy regime scattering is relatively weak and 
perturbation theory converges well. Conversely, in the 
near edge regime (energies less than about 30 eV), the 
range is dominated mostly by the core lifetime 1/T c h 
which is very long for low Z materials or for very shallow 
edges. For those cases an intermediate- or long-range or- 
der theory may be needed, involving multiple-scattering 
to all orders or very large cluster of atoms. 

For polycrystallinc materials or systems with cubic 
symmetry, only the diagonal terms in L and L' survive 
in Eq. (9), corresponding to couplings to various partial 
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local, projected densities of states (LDOS) pi{E), i.e., 

Sr(q,w) « ^(2J + l)|Mi( q)J B)| 2 P i(^)- (14) 

Thus the dynamic structure factor is directly related to 
the LDOS. The coupling terms Mi(q,E) are essentially 
atomic quantities which can be calculated theoretically 
to good accuracy. Thus it may be possible to extract the 
LDOS pi{E) or the density matrix components pl.l'(E) 
from a sequence of experimental measurements of S(q, uj) 
at various q by a suitable inversion procedure of Eq. 
or (|14[) . regarded as a set of linear equations. 

The propagator G = G c + G sc naturally separates into 
intra-atomic contributions from the central atom G c and 
from MS contributions from the environment G sc . Thus 
as in XAS, S(q, u>) can be factored as 



S(q,w)=Sb(9,w)[l + Xq(*)], 



(15) 



where So(q,uj) represents a smoothly varying atomic 
background, and x q (fc) is the "fine structure" due to 
multiple-scattering from the environment. Here we 
have also used a standard XAS notation, where k = 
y/2(u) + Ei) is the photoelectron wave-number. The mo- 
mentum transfer dependence of these two contributions 
to the dynamic structure factor can be analyzed sepa- 
rately as shown in detail below. 



B. Central Atom Contribution 

Much of the qualitative behavior of the spectra can 
be understood in terms of the central atom contribution 
(in the absence of other scatterers), which has been an- 
alyzed for isolated atoms by, e.g., Leapman et al. |3^ . 
This contribution is independent of the direction of the 
momentum transfer, and hence depends only on the mag- 
nitude q. For a condensed system, the states i?i(r, E) = 
i l Ri(r)YL{f), where Ri{r) are radial wave functions and 
Yl{t) spherical harmonics, are scattering states defined 
for the potential of an "embedded" atom in the system. 
However, the deep core states differ little from those of 
isolated atoms and spherical symmetry is still a good 
approximation. This yields the atomic background con- 
tribution 

S (q,u>) = 2(21 + l)\M,(q,E)\ 2 tf(E), (16) 



c 

3 
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FIG. 1: LDOS pf(E) for the central site of fee Al, in the 
absence of other scatterers: I = (solid); I — 1 (long dashes); 
I — 2 (dashes); and I — 3 (short dashes). Note how the 
successive terms p®(E) with increasing angular momentum I 
peak at successively higher energies and overlap each other. 



x r*drRt{r,E)j v (qr)R i {r) 



(17) 



Since ji'(qr) « (qr) near the origin, the dominant term 
for small q is the dipole approximation 1 = 1. The first 
term in a power series expansion of jo(qr) cancels by or- 
thogonality for q — > 0, and hence the first contributing 
terms of jy (qr) from both I' = and I' — 2 behave as 
(qr) 2 . For large momentum transfers such that qa > 1 
and higher excitation energies, successively larger angu- 
lar momentum (I' > 2) excitation channels become more 
important. As a consequence, the behavior of So(q,tv) 
at increasing q reflects that of the LDOS of increasing I. 
This is illustrated in Fig. ^ for the background contribu- 
tion <So(q, ljJ )- F° r EELS however, experimental measure- 
ments of the double differential cross-section d 2 a j dfldui 
still tend to favor the dipole-approximation due to the 
q~ A dependence of the electron scattering cross-section 
(da/dQ) e — e . Thus most of the experimental signal in 
EELS is at small q. 



Fine Structure 



where p® denote the diagonal matrix elements of the den- 
sity matrix. 

The matrix elements can be calculated by rotating the 
system so that the momentum transfer is along z-axis. 
Additionally summing over the m quantum numbers of 
both the initial state and the final state partial wave and 
using properties of 3-j symbols one finds 



|M,(g,£;)| 2 = (2^ + 1)^(2^ + 1) 



U V l 




^From Eq. (9) and (11), the normalized fine structure 
in the dynamic structure factor is given by 



Xq(*) 



Sb(q,w) 



M L (-q, E)p s L ° L ,(E)M L ,(q, E). 



LL' 



(18) 

where pf tL ,(E) = (-l/ir)lmGf tL ,(E). Clearly the direc- 
tional dependence of the spectra on q comes only from the 
scattering contribution, and naturally Xq(^) wm change 
with the magnitude of the momentum transfer due to 
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the q dependence of the coupling terms. This momen- 
tum transfer dependence of the near edge structure has 
been used to study the s ym metries of the core excited 
states in solids H H Q M E IH E H |H |H H| . 
It may also be possible to use this momentum transfer de- 
pendence to obtain information about intermediate- and 
long-range structure, which arc contained in the EXAFS- 
type oscillations X q (fc). 

One way of calculating the scattering contribution to 
the propagator G sc for a big cluster of atoms is to use 
matrix inversion involving the free particle propagator 
G° and the scattering matrix T, i.e., 



G 



L,L' 



,iSi 



[(l-G^T^G^L.L.e 



(19) 



This representation is referred to as full multiple scatter- 
ing (FMS), since it formally includes MS to all orders. 
FMS is often used to treat near-edge spectra where scat- 
tering is strong and the dimensions of the matrix Gl.l' 
are relatively small. At higher energies or whenever the 
MS series converges well, the matrix Gjj & can alterna- 
tively be calculated in terms of a "path expansion," i.e., 
as a sum over all MS paths that a photoelectron can take 
away from the absorbing atom and back |34| . Formally 
the path expansion is given by the sum 



G 



L.L' 



e 



iSi r^iO 



G°TG° + G TG°TG° 



\L,L>e 



iS,, 



(20) 



where the successive matrix products terms represent sin- 
gle, double, and higher order scattering processes. Re- 
markably the path expansion has been found to be gener- 
ally adequate for energies above about 30 eV of threshold, 
where of order 10 2 of the largest amplitude paths suffice 
to yield an accuracy of a few percent for most materials 

m 

Due to the large dimension D s» N{l max + l) 2 of the 
matrix C?lr,i/r', where N is the size of the cluster, exact 
calculations with the path expansion can only be carried 
out for a few, low-order MS paths. To overcome this bot- 
tleneck, an exact separable representation of the propa- 
gator was introduced by Rehr and Albers (RA) |35|. 



G° 



L'R',L"W 



„ikR 

~kR^ 



Y L , iX (kR')Y x , L „(kn"), (21; 



where R = |R' — R"| is the bond distance. The gen- 
eralized spherical harmonic expansion coefficients Yl',\ 
and 1a, L" converge rapidly in powers of 1/kR, which is 
always <C unity for unoccupied states fc > hp, and hence 
this representation can be severely truncated. For FMS 
calculations the RA representation of the propagators in 
Eq. (|21() is stable and converges rapidly. For the path ex- 
pansion, the RA approach is exact for single scattering, 
and typically only six terms in A suffice to give accu- 
racies to of order a percent or better over the typical 
range of wave numbers encountered in XAS experiment, 
k F < k < 20 A" 1 , i.e., E F <E < 1500 cV. 

With the separable representation, one can sum over 
all intermediate angular momentum variables (l,m) at 



each site and define local scattering matrices. Then the 
contribution G°TG°T ■ ■ ■ G° to the total propagator from 
a given A-atom path T = [Ri, R2, . . . Rjv = R-o] can be 
factored as a matrix product over small (typically 6x6) 
matrices Hal 



L 



which is the analog of the scattering amplitude f(9) in 
plane-wave scattering theory, and a termination matrix 
™£,L' N (Pi, Pn) = Ylm(Pi)Yl n ,\ n (p n ), where p = fcR 
[35|. For large I or L', it may be necessary to increase 
the dimension of A in the termination matrices m 



Ai ,Aj\ 
L.L' 



beyond 6 for adequate con verg ence. Specifically we ob- 
tain for each path (cf. Ref. |3q1. 



G T 



D i(pi+P2 + ---PN+Sl+S l l) 



L,L' 



E m ^* '(Pi'Pn) X 

{A} 



P1P2 ■■■ PN 

(PiV>Piv-i)---^A 2j A 1 (p 2J P 1 )] -(22) 

In our code, the dependence on the bond vectors p and 
p' is simplified using rotation matrices and Euler angles. 
This expression is the similar to one used for analysis of 
XAS using path expansion Q] . In contrast to the case for 
XAS which is dominated by the dipole-approximation, 
however, all angular momentum channels now contribute 
with q dependent couplings. Thus the fine structure from 
the path expansion can be written as 



Xq(*) 



where 



-Im^ 



i(pi+p 2 + ---pjv) 



PlP2 ■ ■■ PN 



{A} 



'(pi.Piv>q) 



[Fx n ,x n _ 1 (p n ,Pn-i) ■ ' • i? A 2 ,A 1 (p 2l Pi)] ■ (23) 



# Al ' A ~(p iPjv ,q) 



E 

L,L> 



M L (-q)m x L ^ N ( Pl ,p N )M LI (ci). 



(24) 



Note that only the coupling terms H Xi ' Xn (p 1 ,p N , q) de- 
pend on q, while the scattering contribution is a product 
of low order scattering matrices F\.A' which is dependent 
only on the material. When the dominant contribution 
comes from the dipole limit I = I' = 1, the dependence on 
q has the form (q • Ri)(q - Rat). This form emphasizes 
paths beginning or ending in the direction q, and acts 
like a "search-light" in probing the structure of a sys- 
tem. For polycrystallinc materials or for measurements 
averaged over all q, this dependence averages out. Simi- 
larly, the contributions from the higher angular momenta 
couplings probe other symmetries. 

By defining a q-dependent effective scattering ampli- 
tude jeff (q, k) one can use the path expansion (Eq. (12 3|) ) 
to express the fine structure Xq(^) m a form identical to 
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the original EXAFS equation, 

*»(*) ^oE lfc l% k)l M2kR+^ k )e-^e-^ k \ 

r (25) 

where k = y/2(uj + E t ), R = (1/2)2^ is the ef- 
fective path length, and the prefactor Sq « 0.9 is a 
many-body amplitude factor which accounts for inelas- 
tic losses (satellite excitations) beyond the quasi-particle 
approximation!^. Here the effects of thermal and struc- 
tural disorder are included using a configurational aver- 
age of Eq. These damping effects can be approxi- 
mated by a Debye- Waller factor exp(— 2a 2 k 2 ), where a 2 
is the correlated mean-square variation ((SR) 2 ) for each 
scattering path pj, while anharmonic terms from the 1st 
and 3rd cumulants are generally weaker and contribute 
to the phase. Typically a 2 is of order 10~ 2 -10 -3 A 2 . 
For the case of full multiple-scattering, these effects can 
be calculated to a good approximation by including a 
similar Debye-Waller factor exp(— a 2 k 2 ) in each propa- 
gator G lr,l'R' (k) ■ Close to an absorption edge where 
k is small, the Debye-Waller factors are of order unity 
and can be ignored. In practical calculations cr 2 (T) can 
be approximated to reasonable accuracy using the corre- 
lated Debye model evaluated for the Debye temperature 
Q]j of a material pj. 

III. EXAMPLES 

A. Be A'-edge 

As mentioned above the momentum transfer magni- 
tude q and directional q dependence of the K-edge in 
Be metal was measured and analyzed by Sternemann et 
al. [3. They found that the main changes in the spec- 
tra with increasing q could be explained by the increas- 
ing contribution from monopole (s-to-s) type excitations. 
This leads to a net decreasing anisotropy, i.e. decreasing 
dependence on q of the edge, since s-type states natu- 
rally have no directional dependence. This conclusion is 
supported by ground state local density of states calcu- 
lations using full-potential linear-augmented plane-wave 
and excited state calculations using the BSE method of 
Ref. 

In Fig. |2 we compare the inelastic scattering spectra 
calculated with the approach of this paper for the Be 
AT-edge with the experimental results of Ref. ■ From 
this comparison it is clear that our present calculation is 
capable of reproducing both the momentum transfer di- 
rection and magnitude dependence of the dynamic struc- 
ture factor. However, there are noticeable discrepancies 
in some fine details in the spectra within the first 15 eV. 
For the momentum transfer perpendicular to the c-axis, 
the 120 eV features are more pronounced in the RSMS 
calculation than in the experiment for both the low- and 
high momentum-transfer. However, the main differences 
between these spectra is well reproduced by present cal- 
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FIG. 2: Comparison between experiment [iH (top panels) 
and inelastic loss spectra calculated using the BSE (middle 
panels) and the RSMS method of this work (lower panels) for 
the Be A'-edge. The spectra are shown for q perpendicular to 
the c-axis (left panels), and q along the c-axis (right panels). 
Spectra are shown for two values of momentum transfer q as 
indicated in the figure labels. The theoretical spectra was 
shifted by 4.5 eV to align with the experimental edge. 



culation. In the middle panel of Fig. [2 we show for com- 
parison the results from the BSE method |18l l29|. Con- 
sidering how sensitive the problem of calculating core- 
excitation is, the agreement between the two different 
theoretical methods is rather good. This serves to val- 
idate the one-electron, final state rule calculations, at 
least for this case. As noted above, most of the differ- 
ences occur rather close to the edge, which is also the re- 
gion most sensitive to the details of the core hole-electron 
interaction and to the nature of the scattering potentials. 
Beyond the first 15 eV the agreement between the two 
methods is very good. For q along the c-axis, agree- 
ment between the two calculations is also rather good 
for the energy range shown. For q perpendicular to the 
c-axis, however, the comparison is somewhat mixed. For 
the first peak in the spectra it appears that RSMS accu- 
rately predicts the experimental result. The difference is 
that the first peak becomes a shoulder at low momentum 
transfer and a peak at higher momentum transfer. On 
the other hand the details of the spectra between 120- 
130 eV appear to be slightly better reproduced by the 
BSE. 

As noted in Sec. II. C, the momentum transfer depen- 
dent changes in the spectra can be understood in terms 
of the changing weights of the different components of 
Pl.l'(E) in the single crystal case, and the LDOS pi(E) 
for polycrystalline systems. In Fig. [21 we show the s- p- 
and d-DOS for the core excited state at the site of the core 
hole in Be metal. Also shown are the calculated inelastic 
scattering spectra for two values of momentum transfer 
along the x-axis. The calculated spectra are shifted so 
that the Fermi-energy is aligned with that obtained in 
the LDOS calculation, where E = is the vacuum level. 
On this scale the Fermi energy for Be as calculated by our 
RSMS code is E F = -9.4 eV. The p-type LDOS close to 
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FIG. 3: Local projected DOS for Be together with inelas- 
tic scattering spectra for two values of momentum transfer q 
along the a:-axis. The Fermi level is at -9.4 eV. 



the Fermi energy is composed of two peaks at —8 eV and 
+4 eV and a broad shoulder at —4 eV. In the small mo- 
mentum transfer regime the dipole limit is reached, and 
the general shape of the p-type LDOS clearly shown in 
the shape of the calculated spectrum. The s-type LDOS 
on the other hand is composed of three sharp peaks at 
—8, —4 and +3 eV. Compared to the low q spectrum, the 
high momentum transfer spectrum has sharp peaks at 
these same energies. The coincidence of the peaks is due 
to hybridization between the various angular momentum 
components. Clearly the changes in the near edge struc- 
ture with increasing q can be attributed partly to the 
increasing contribution from the s-to-s type transitions. 
Also the changes due to the matrix elements that weight 
the different excitations have to be considered when com- 
paring LDOS to experimental spectra. This can clearly 
be seen when comparing low momentum transfer spectra 
along either the x- or c-axes. Although both exhibit a p- 
type LDOS structure, the weights of the different features 
are strongly dependent on the magnitude and direction 
of the momentum transfer. Also the overlap of the dif- 
ferent types (s vs p) of LDOS is clearly demonstrated in 
these results. The relative contributions to the sum of 
Eq. (|14f) from monopole (s-to-s), dipole (s-to-p) transi- 
tions and quadrupole (s-to-d) are illustrated in Fig.0]for 
the higher momentum transfer q = 4.78 a.u. along the c- 
axis. At lower momentum transfer the contributions from 
the monopole transitions and quadrupole transitions are 
negligible. At this relatively high momentum transfer 
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FIG. 4: Comparison of contributions to the inelastic loss 
spectra from monopole and dipole transitions, for momentum 
transfer q — 4.78 a.u. along the c-axis. 

the dipole forbidden s-to-s transition is comparable in 
magnitude to the dipole allowed s-to-p transition. The 
s-to-d contribution is small at low energies but becomes 
increasingly important as the energy is increased, again 
reflecting the energy dependent weight of different exci- 
tation channels. 

The total spectrum changes also as function of the 
magnitude of the momentum transfer because of the 
changes in the channel specific spectra. This change is 
due to the excitation energy E dependence of the final 
state partial wave Rl(?, E). As the momentum transfer 
changes the matrix elements Mr,(q, E) = (i?i(i?)|e lq ' r |i) 
make the channel (i.e., L) specific spectra change at 
an energy E dependent rate. In Fig. [S] the momentum 
transfer dependence of the dipole allowed and monopole 
(dipole forbidden) transitions are examined separately. 
The spectra are scaled so that the first peak for both 
momentum transfers is the same height. Besides this en- 
ergy independent scaling factor there is also a momentum 
transfer dependence on the shape of the channel specific 
spectra. At a qualitative level this change is mostly vis- 
ible in the central atom contribution. The momentum 
transfer dependent changes in the fine structure are less 
noticeable. 



B. Al Li edge 

As a final example we examine the core excited states 
in fee Al. The calculated x-ray absorption cross-section 
and the dynamic structure factor for Li-edge at different 
momentum transfers are compared in Fig. To make 
contact with S*(q, u>) the absorption cross-section a a b s {u) 
was divided by the excitation energy, since in the limit 
q^0, i 

dabs (w) OC U)S(q,LO). (26) 
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FIG. 5: Changes in the shape of different contributions to 
the Be JT-edge spectra, as a function of the magnitude of the 
momentum transfer along the c-axis for two values of q listed 
in the figure labels. The upper panel shows the changes for 
the monopole contribution and the lower panel, those for the 
dipole allowed transitions. The spectra were normalized so 
that the height of the first peak is the same for both values 
of q. 
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FIG. 6: Comparison of inelastic scattering and x-ray absorp- 
tion spectra at the Al Li-edge using the path expansion for 
several values of momentum transfer as indicated in the figure 
labels; the curve a(uj)/uj is the result for q — 0. For clarity, 
the inelastic scattering spectra was slightly shifted vertically. 
All the spectra where scaled to be shown on the same figure. 
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FIG. 7: Fine structure Xq(^) f° r Al at the Li edge for different 
values of momentum transfer as listed in the figure labels, 
together with the x-ray absorption fine structure x(k), i.e., 
q = 0. The curves are shifted for clarity. 



Comparison of a a b s {uS)/u) and S{t\,w) for q = 0.24 a.u. 
clearly shows the well known and experimentally veri- 
fied [4| fact that inelastic scattering results at low mo- 
mentum transfer can be used to obtain the same EX- 
AFS information as in XAS. Since we are interested in 
discussing the general trends of the momentum transfer 
dependence of Si(q, oS) in the extended energy range we 
have not included the Al L 2 ,3-edge (edge located at about 
80 eV) to this theoretical demonstration. For comparison 
with momentum transfer dependent experiments for this 
energy range (we are not aware of such experiments in 
any material) one would need to include these edges. As 
q is increased to q = 1.5 and 4.5 a.u., the overall shape of 
the spectra changes. As mentioned above, this change is 
mostly due to changes in the central atom contribution. 
Although the EXAFS-likc oscillations appear to dimin- 
ish with higher momentum transfer, most of the broad 
oscillatory structure remains. 

In Fig. {7\ we examine the momentum transfer depen- 
dence of the q-dependent fine structure Xq(fc)- Here we 
show the XAS fine structure x(fc), together with x q (&;) 
for several values of q at the Al Li-edge. Again it is clear 
that for low q the XAS result is reproduced. When the 
momentum transfer is increased, however, the shape of 
the fine structure of the spectra is modified in a way that 
cannot be explained by simple scaling factors, due to the 
mixture between the various angular momentum contri- 
butions. The changes are perhaps strongest close to the 
edge (at small k), but even the high k spectra is modi- 
fied. This suggests that it may be desirable to decompose 
the spectra into the various angular momentum com- 
ponents Xz(/c) prior to additional analysis. This is fur- 
ther illustrated in Fig. [5] which shows the Fourier trans- 
form of the EXAFS Xq m -R-spacc, phase corrected by 
the dominant central atom p-wave phase shift exp(2i(5i). 
Such Fourier transforms have peaks close to the near- 
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FIG. 8: Fourier Transform of Al fine structure for several 
values of momentum transfer as listed in the figure labels. The 
differences reflect the interference between different channels 
contributing to the spectra. 

neighbor distances and provide a characterization of the 
near neighbor radial structure. Thus they can be used 
for quantitative fits of EXAFS in i?-space. Note that for 
small q the transform is insensitive to q but becomes sub- 
stantially more complex for the larger values due to the 
overlapping contributions from the various channels. 

IV. CONCLUSIONS 

We have shown how the real-space multiple-scattering 
(RSMS) approach can be applied to model non-resonant 
inelastic scattering from deep core levels for arbitrary 
condensed systems over a broad spectral range. The 
approach is a generalization to finite momentum trans- 
fer of that used to model XAS based on the indepen- 
dent electron approximation and the final state rule. In 
contrast to XAS where dipole selection rules apply, cou- 
plings to all angular momentum components are impor- 
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